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ABSTRACT 

Context. The complexity of the common envelope phase and of magnetic stellar wind braking currently limits our 
understanding of close binary evolution. Because of their intrinsically simple structure, observational population studies 
of white dwarf plus main sequence (WDMS) binaries hold the potential to test theoretical models and constrain their 
parameters. 

Aims. The Sloan Digital Sky Survey (SDSS) has provided a large and homogeneously selected sample of WDMS binaries, 
which we are characterising in terms of orbital and stellar parameters. 

Methods. We have obtained radial velocity information for 385 WDMS binaries from follow-up spectroscopy, and for 
an additional 861 systems from the SDSS sub-spectra. Radial velocity variations identify 191 of these WDMS binaries 
as post common envelope binaries (PCEBs). Orbital periods of 58 PCEBs were subsequently measured, predominantly 
from time-resolved spectroscopy, bringing the total number of SDSS PCEBs with orbital parameters to 79. Observational 
biases inherent to this PCEB sample were evaluated through extensive Monte Carlo simulations. 

Results. We find that 21 — 24% of all SDSS WDMS binaries have undergone common envelope evolution, which is 
in good agreement with published binary population models and high-resolution HST imaging of WDMS binaries 
unresolved from the ground. The bias corrected orbital period distribution of PCEBs ranges from 1.9 h to 4.3 d and 
follows approximately a normal distribution in log(P or b), peaking at ~ 10.3 h. There is no observational evidence for a 
significant population of PCEBs with periods in the range of days to weeks. 

Conclusions. The large and homogeneous sample of SDSS WDMS binaries provides the means to test fundamental 
predictions of binary population models, and hence to observationally constrain the evolution of all close compact 
binaries. 

Key words, stars: novae, cataclysmic variables - stars: binaries: close - stars: white dwarf-main sequence binaries - 
stars: pre-cataclysmic variables - stars: post common envelope binaries 



1. Introduction 

It is well known that all types of close compact binaries, 
including X-ray binaries, cataclysmic variables, and dou- 
ble degenerates, form through common envelope evolution 



(jPacz vhski 1976). This phase is required to explain the ob- 
served small binary separations in compact binaries. The 
main concept is that, for a suitable range of orbital separa- 
tions, the massive star (i.e. the primary) will fill its Roche 
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lobe as it evolves up the giant (or asymptotic giant) branch, 
and will then in most cases become susceptible to dynami- 
cally unstable mass transfer. The transferred material can- 
not cool as fast as it is being transferred and the core of 
the primary star and the companion star (also referred to 
as secondary star) are both engulfed by the envelope of 
the former. The envelope is not co-rotating with the bi- 
nary, creating drag forces which transport orbital angular 
momentum and energy from the binary to the common en- 
velop e, leading to a shrinkage of the orbit (jRicker fc Taaml 
120081) . Eventually, the common envelope is ejected, leaving 
behind a post common envelope binary (PCEB). 

Although the basic concept of the common envelope 
phase has been outlined 30 years ago, it is still the most 
poorly understood phase of compact binary evolution. 
While we h ave learnt that this phase must be very short 
(< 10 3 yrs) (|Hiellming fc Taamlll99ll: IWebbinkll2008jb cur- 
rent hydrodynamical simulations of common envelope evo- 
lution are still unable to unambiguously link the initial pa- 
rameters of the binary with the outcome of the common 
envelope, and numerically too expe nsive to be run on larg e 
ranges of initial binary parameters (Packer & Taam 2008) . 
Therefore theoretical binary population synthesis studies 
(BPS) usually adopt a simple parametrisation: a certain 
fraction («ce) of the binary's binding energy which is re- 
leased in the spi ralling- in proce ss is used to unbind the com- 
mon e nvelo pe (fWebbinkl Il984f ) . How ever, iNelemans et al.1 
(|2000D and INelemans fc Tout! (|2005D find difficulties even 
with this scaled energy relation when trying to reproduce 
the evolution of double white dwarf binaries, and propose 
a prescription based on angular momentum conservation. 
A fundamental consequence of this approach is that clear 
observational constraints on these parameters are necessary 
to make any meaningful predictions on the Galactic popu- 
lation of compact binaries. 

PCEB^] composed of a white dwarf and a main- 
sequence (WDMS) companion represent the most promis- 
ing population to derive such observational constraints, 
as they are the intrinsically most abundant population of 
PCEBs and hence easily accessible to intense studies using 
the current suite of 2-8 m telescopes, they underwent only 
one common envelope phase, are not altered by ongoing ac- 
cretion, and their stellar components are structurally sim- 
ple. In general, WDMS binaries comprise two fundamen- 
tally different types of systems: wide binaries that evolved 
like single stars, i.e. without interaction, and PCEBs. The 
major shortcoming of using PCEBs for studying compact 
binary evolution h as been their small number. The com- 
prehensive study of lSchreiber fc Gansickd ([20031 ) contained 
only 30 systems with accurate binary parameters, and high- 
lighted that the sample of PCEBs known at that time was 
not only small but also heavily biased towards hot white 
dwarfs and late secondary star spectra l type s. This bias 
was identified by ISchreiber fc Gansickd (2003) as a natu- 
ral consequence of the fact that most of these PCEBs were 
discovered in blue-excess surveys. 

Over the p ast few years, t he Sloan Digita l Sky 
Survey (SP S S. lYork et al.1 [20001 lAbazaiian etafl [20091: 
lYannv et al.l I2009T) has totally changed this situation. 
Thanks to its broad ugriz colour space, and mas- 
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sive spectroscopic foll ow-up, SPSS has i dentified over 
1900 WPMS binaries (tSmolcic et al l 12004 ISilvestri et al 
1 20061 I2007t ISchreiber et alJ l2007t iHeller et al.l 1200' 
iRebassa-Mansereas et al.l l2007. 201Q). We have initiated an 
intensive follow-up study of the SPSS WPMS binaries, with 
the ultimate aim to substantially improve our understand- 
ing of compact binary evolution. First results on individ- 
ual system s and small sub-sets of this sample have been 



1 For the remainder of the paper, we will use PCEB as a 
synonym for a short-period WDMS binary that underwent a 
common envelope phase. 



published (ISchreiber et ail 20081: IRebassa-Mansergas et all 
1 20081: iNebot Gomez-Moran et al.ll2009l: iPvrzas et all 12009: 
ISchwope et al.l |2009() . as well as several studies that 
analyse the t otal sample in the context of close bi- 
nary evolution (ISchreiber et aifeOlOHZorotovic et al.ll2010l: 
IRebassa-Mansergas et al.ll2011f) . Here we present one of the 
fundamental results of this project: the orbital period dis- 
tribution of 79 SPSS WPMS binaries, which we discuss in 
the context of current theories of compact binary evolution. 



2. Overall strategy and input sample 

Our aim of characterising a large and homogeneous sample 
of PCEBs requires a multi-stage approach. The first step 
is to identify all WPMS binaries within the SPSS spec- 
troscopic data release. The latest catalogues of WPMS bi- 
naries from SPSS list 190 3 WPMS systems, - 1600 from 
IRebassa-Mansergas et al.1 (|2010T > based on the SPSS Pata 
Release (PR) 6, plus another ~ 30 that have been dis- 
covered by us ([Schreiber et al.l [200711 as part of the Sloan 
Extension for Galactic Und erstanding and Exploration 
(SEGUE, lYannv et all I2009D . These catalogues include 
white dwarf temperatures and masses, as well as the spec- 
tral types of the companion stars. In brief, stellar parame- 
ters are derived from fitting the SPSS spectrum to template 
spectra made of white dwarf plus dM templates. The spec- 
tral type of the secondary star is derived from the best fit 
and the mass is calculated using a spectr al type-mass rela- 
tion (see IRebassa-Mansergas et al.1 12007L for details) . The 
mass and the temperature of the white dwarf are derived 
from fitting both lines and continuum (after the contribu- 
tion from th e companion has be en subtracted) to stellar 
models from iKoester et aD (|2005h . 

The second step is to identify the PCEBs among the 
entire sample of W PMS binaries. The com prehensive pop- 
ulation synthesis bv lWillems fc KoTbl (|2004l ) suggested that 
~ 75% of all Galactic WPMS binaries are wide systems, 
and that the vast majority of the remaining < 25% un- 
derwent a common envelope phase, hence several hundred 
PCEBs are to be expected among the SPSS WPMS sam- 
ple. A small fraction of WPMS binaries appear as blended 
white dwarf plus M-dwarf pairs in the SPSS imaging, and 
were fl agged as wide systems by IRebassa-Mansergas et ail 
([2010D . However, even those WPMS binaries that are gen- 
uinely unresolved on the SPSS images can still have orbital 
separations of many tens of astronomical units. The key 
to identify the close WPMS binaries is hence to carry out 
a radial velocity survey of the spatially unresolved SPSS 
WPMS binaries, obtaining at least two spectra, separated 
by at least one night, of as many WPMS binaries as pos- 
sible. In a final step, intense follow-up spectroscopy of the 
PCEBs, i.e. WPMS binaries found to be radial velocity 
variable, is needed to measure their orbital periods. 
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3. Observations, reduction, and radial velocity 
measurements 

We obtained spectroscopic follow-up observations for 385 
systems (Sect J3.ip . and used exclusively SDSS spectro- 
scopic data for another 861 systems, bringing the number 
of SDSS WDMS binaries with at least two reliable radial 
velocity measurements to 1246. In addition, we obtained 
photometric time-series for seven systems (Sect l3.3j) . A log 
of the observations is given in the online version of this 
article (Table [TJ, listing the name of the object, the ugriz 
magnitudes, the telescopes, the number of spectra taken 
and the duration of the observations. 



3.1. Spectroscopy 

Throughout the period January 2006 to July 2010, we have 
used the 3.5 m telescope at Calar Alto (CA), the 3.5 New 
Technology Telescope (NTT), the 4.2m William Herschel 
Telescope (WHT), the two 6.5 m Magellan telescopes 
(Baade/Clay), the two 8m Gemini telescopes (GN/GS), 
and the 8 m Very Large Telescope (VLT) to obtain spec- 
troscopy of 385 SDSS WDMS binaries. Details on the in- 
strumentation and setup are given in Table [5] 

The data obtained at Calar Alto were reduced within 
MIDAS0 a nd the spect ra were extracted using the optimal 
algorithm (|Horndll986tl . The data from all other telescopes 
was reduced using the STARLINK0 packages FIGARO and 
KAPPA, spe ctra were optimally extracted using PAMELA 
(|Marshlll989l ) and wavelength-calibrated within MOLLIS- 
For all the spectra, the wavelength calibration was checked, 
and if necessary corrected, for telescope/instrument flexure 
using the night sky lines. The spectra were flux calibrated 
and corrected for telluric absorption using observations 
of the stars BD+28°4211, BD+33+°2642, BD+25°4655, 
Feige66, FeigellO, GD108, and LTT3218. For additional 
details on the r educti on and calibration of the d a ta, se e 
Schreiber et al.l (120081); IRebassa-Mansergas et all (|2008jk 
Nebot Gomez-Moran et all (|2009t) . 

In addition to our own follow-up spectroscopy, we 
made use of the spectroscopic data from SDSS. Every 
SDSS spectrum is the average of (typically) three 15min 
exposures. These sub-spectra are publically available 



Table 2. Observational technical details. 



pt 

since DR6 (lAdelman-McCarthy et all [20081). As shown b; 
IRebassa-Mansergas et al.1 ( 20081 ) and lSchwope et al.l (1 200' 
it is possible to use many of these sub-spectra to mea- 
sure radial velocities and identify short-period WDMS bi- 
naries. We retrieved the subspectra for all 1903 systems 
in our input catalogue. Discarding spectra of insufficient 
signal-to-noise ratio, we were able to measure at least one 
radial velocity for 1147 WDMS binaries, and we measured 
a total of 5171 radial velocities using the SDSS data. The 
wavelength calibration of the SDSS spectra can be accu- 
rate to ~ 2kms~ 1 in dark nights, while in bright nights 
there can be shifts up to 20 k ms" 1 . A systemat i c erro r of 7 
kms" 1 has b een re ported bv IStoughton et all ((2002) and 
lYannv et all (|2009h . Although the SEGUE spectroscopic 
data has been corrected for that offset, this shift has only 
been applied to the combined spectra and not the indi- 
vidual subspectra. We take this into account by adding in 
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opectroscopy 








Tel. /Inst. 


Slit 


Spcct . 


O \J ee L . 




it 01 


width[' 


range [A] 


r s i 

res. [A 


[km s J 


RVs 


VLT/FORS2 


1 


7830-9570 


2.2 


7.1 


334 


Gcmini-N/GMOS 


1 


7770-8470 


3.5 


7.1 


124 


Gomini-S/GMOS 


1 


7770-8470 


4.1 


8.6 


26 


Clay:f/LDSS3 


0.75 


5800-9980 


4.8 


10.7 


302 


BaadcJ/IMACS 


0.75 


7700-8560 


1.5 


11.0 


187 


WHT/ISIS 


1 


7600-9000 


1.6 


9.2 


524 


NTT/EMMI 


1 


7770-8830 


2.8 


11.4 


284 


NTT/EFOSC2 


1 


7000-8600 


6 5 


19.5 


94 


CA3.5/TWIN 


1.5 


7500-8500 


1.6 


11.6 


443 


CA3.5/MOSCA 


1.5 


5400-9000 


7 


17.0 


12 


SDSS 


3+ 


3900-9100 


4 


19.1 


5168 






Photometry 








Tel. /Inst. 




Resolution ["/pix. 




# of images 


IAC80/CAMELOT 




0.6 




397 




DuPont/CCD 




0.259 




337 




CA2.2/CAFOS 




0.53 




1083 





Notes. 1 " SDSS uses a fibre-fed spectrograph. * The two Magellan 
6.5 m telescopes. 



quadrature 10 kms 1 to the statistical errors on the radial 
velocities measured from Gaussian fits (see Sect. 13. 2|) . 

3.2. Radial velocity measurements 

We measured the radial velocities of the companion stars 
by fitting a second-order polynomial plus two Gaussians 
to the observed 8183/8194 A Nal absorption doublet. The 
separation of the two Gaussians was fixed to the labora- 
tory value of the Nal doublet. We also measured the ra- 
dial velocities from the Ha emission line (where present) 
by fitting a second-order polynomial to the underlying 
continuum, plus a Gaussian for th e emission line (see 
IRebassa-Mansergas et al.|[2007l 120081 for further details on 
the method). Radial velocities, errors, time of the obser- 
vations and the telescope used for the identification of the 
system are available via the CDS; an excerpt is presented 
in Tabled 

3.3. Photometric follow-up observations 

We carried out time-series photometry for 7 WDMS bina- 
ries to probe for orbital variability (Tablc[TJ where the avail- 
able spectroscopy suggested a short orbital period. These 
data were obtained between September 2006 and August 
2009 using the 2.5 m DuPont telescope, the 2.2 m telescope 
at Calar Alto, and the 0.8 m IAC80 telescope. All obser- 
vations were carried out with red filters (Bessel R, Sloan r 
or i) to maximise the sensitivity to changes in the bright- 
ness of the secondary star (ellipsoidal modulation and/or 
reflection effect). The IAC80 data were reduced using stan- 
dard packages in IRAI-E The DuPont and Calar Alto data 
were de-biased and flat-fielded in a standard fashion within 
MIDAS, an d instrumental magnitu des were derived using 
SExtractor (jBertin fc Arnoutsl ll996') . For a full description 
of the photometric reduction pipeline see iGansicke et al.l 
(2004|). 



2 http:/ /www. eso.org/sci/software/esomidas 

3 starlink.jach.hawaii.edu/starlink 



www 



warwick.ac.uk /go /trmarsh / software 



5 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under contract with the 



National Science Foundation, http://iraf.noao.edu 
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Fig. 1. Upper panel: frequency of the number of spec- 
tra taken per system from our spectroscopic campaign 
(dashed), from the SDSS (solid) and total combined num- 
ber of spectra (gray). Bottom panel: distribution of time 
span between the first and the last spectrum taken for iden- 
tification of close and wide binaries. See text for a detailed 
explanation. 
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Fig. 2. Distribution of the probability of measuring strong 
radial velocity variations among the WDMS binaries in the 
sample. Those above a 3a (P(x 2 ) > 99.73) threshold are 
considered to be PCEBs. Line coding as in Fig. [TJ A zoom 
in the region between 0.995 and 1 is shown in the box, note 
that the bin size is smaller. For a detailed explanation see 
text. 



The results of the photometric observations are dis- 
cussed in Sect. 14. 2. T1 and 14.2.31 and details on individual 
systems are given in the online version ( Appendix lA"]) . 

4. Observational results 

4.1. Post common envelope binary identification 

Our observational program started back in 2005. The orig- 
inal strategy for identification of PCEBs was to take two 
to three spectra for each WDMS binary at random times, 
with the only constraint that they should be separated by 
at least one night. Any system displaying significant radial 
velocity variations was then flagged for intense follow-up 
spectroscopy to measure its orbital period. 

With the wealth of SDSS increasing with every 
data release, our strategy ha s evolved. Within DR5, 
iRebassa-Mansergas et all (|2007l) found that - 100 WDMS 
had multiple SDSS spectra, and identified a certain frac- 
tion of PCEBs among them. As mentioned in Sect. l3.1l each 
SDSS spectrum is typically the average of three individual 
exposures. Occasionally, the subspectra were obtained on 
different nights, and for a small number of objects more 
than three subspectra are available. These subspectra were 
made publically available from DR6 onwards, and we made 
use of them for the identification of PCEBs. For ~ 100 sys- 
tems, the signal-to- noise ratio of the subspectra was insuffi- 
cient for an accurate radial velocity measurement. In these 
cases, we combined the radial velocity determined from the 
co-added SDSS spectra with our own measurements. For 
those WDMS binaries that have SDSS data obtained on 
different nights, we usually took only one additional spec- 
trum to have an independent measurement. 

Figured shows the distribution of the number of spectra 
per system used to separate PCEBs and wide binary can- 
didates (upper panel), and the time elapsed between the 



first and last spectrum for a given object (bottom panel), 
with the SDSS data (solid line), our follow-up data (dashed 
line), and the combined data (gray). Most of the WDMS 
binaries have 3 SDSS subspectra which were taken in the 
same night, and are hence sensitive only to short orbital 
periods. A second peak around a few days indicates SDSS 
repeat observations of the same spectroscopic plate. Our 
own spectroscopy was typically carried out with a separa- 
tion of ~ 1 — 3 days in visitor mode observing runs, with 
a long tail coming from service mode observations on the 
VLT and on Gemini. We make the reader note that the 
combined distributions are not the direct sum of the other 
two, there are systems that were observed only once by us 
but had several SDSS subspectra, or the other way round. 

For the identification of radial-velocity variable stars, 
i. e. PCEBs, we analysed the radial velocities measured 
from the Na I doublet and the Ha emission line separately, 
since these two velocities can differ by more than the statis- 
tical errors (jRebassa-Mansereas et al.ll2007l ). e.g. if Ha orig- 
inates predominantly in the inner he misphere of the com- 
panion. We followed the approach of iMaxted et al.l ((2000) 
to evaluate the significance of the measured radial velocity 
differences. We calculated for each system the x 2 statistic 
with respect to the mean radial velocity, which represents 
the best fit for a non-variable system. For the hypothe- 
sis of a constant radial velocity, the probability Q of ob- 
taining a given \ 2 is high when radial velocity variations 
are dominated by random noise, and low for intrinsically 
radial-velocity variable systems. The probability of a sys- 
tem showing high radial velocity variations is hence given 
by P(x 2 ) = 1 — Q. We consider strong PCEB candidates 
those systems showing radial velocity variations at a > 3er 
significance, or P{x 2 ) > 99.73%. The number of spectra 
used for each system, time span covered by these spectra, 
and the probabilities P(x 2 ) are available via the CDS; for 
an extract see Table |U 
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Table 3. 

binaries. 



Radial velocities used for the identification of the 



Name 


RV 


err 


hjd 


Teles. 


SDSSJ 


[kms -1 ] 


[kms~ r 


[days! 




002749.99 


- 001023.3 


13.2 


15.3 


2452262.628080 


SDSS 


002749.99 


- 001023.3 


11.0 


17.3 


2452262.640324 


SDSS 


002749.99 


- 001023.3 


7.0 


17.0 


2452262.652580 


SDSS 


002749.99 


- 001023.3 


8,1 


10.5 


2454347.439589 


WHT 


002749.99 


- 001023.3 


15.8 


6.7 


2454348.551886 


WHT 


002801.68 


+ 002137.7 


-25.6 


22.4 


2452262.628105 


SDSS 


002801.68 


+ 002137.7 


1.8 


21.0 


2452262.640349 


SDSS 


002801.68 


+ 002137.7 


0.0 


33.6 


2452262.652605 


SDSS 


002801.68 


+ 002137.7 


-4.8 


5.0 


2454381.694818 


VLT 


002801.68 


+ 002137.7 


1 Q 

-o.o 


. u 


2454382.804165 


VLT 


002950.16 


+ 003225.8 


18.3 


21.0 


2452262.628153 


SDSS 


002950.16 


+ 003225.8 


6.2 


17.8 


2452262.640397 


SDSS 


002950.16 


+ 003225.8 


31.5 


27.4 


2452262.652652 


SDSS 


002921.27 


- 002032.8 


-62.6 


28.0 


2452262.628107 


SDSS 


002921.27 


- 002032.8 


-25.6 


17.0 


2452262.640351 


SDSS 


002921.27 


- 002032.8 


-22.3 


18.7 


2452262.652607 


SDSS 


003221.86 


+ 073934.4 


-326.6 


8.3 


2454681.908960 Gcmini-S 


003221.86 


+ 073934.4 


-27.9 


5.9 


2454687.871683 Gcmini-S 


003336.49 


+ 004151.3 


67.8 


17.2 


2452930.850863 


SDSS 


003336.49 


+ 004151.3 


90.9 


16.4 


2452932.782491 


SDSS 



It is clear that any radial velocity survey will be subject 
to observational biases related to the spectral resolution 
employed, the number of spectra obtained for each object, 
and the temporal spacing of these spectra. These biases 
are evaluated in Sect. 15.21 They are also responsible for the 
differences between the probability distribution from the 
SDSS data and from our own follow-up observations, which 
are discussed in Sect. [5] 



4.2, Orbital period measurements 

As mentioned above, 191 WDMS binaries exhibited 
significant radial velocity variations and were flagged 
for detailed follow-up observatio ns. 15 orbit al per i- 
ods were alre a dy p u blished by iPvrzas et al.l (12011 ); 
ISchreiber et all (120081 ); iRebassa-Mansergas et all (12008); 



iNebot Gomez-Moran et al.l (|2009) and IPvrzas et al.l 
(|2009D . Here we report period determinations for 58 addi- 
tional systems, four of which were measured from follow-up 
photometry, the rest from follow-up spectroscopy. 



Table 4. Identification statistics. 



Name 


N 


N obs 


At [days] 


Prob. 


SDSSJ 000504.91 +243409.6 


10 


3 


294.821 


0.74189 


SDSSJ 000531.09 - 054343.2 


6 


1 


427.720 


1.00000 


SDSSJ 000935.50 + 243251.2 


16 


9 


377.808 


0.97822 


SDSSJ 001749.24 - 000955.3 


2 


2 


3.935 


1.00000 


SDSSJ 002620.41 + 144409.5 


3 





0.029 


0.32658 


SDSSJ 003221.86 + 073934.4 


2 


2 


5.963 


1.00000 


SDSSJ 004020.07 + 154156.3 


4 





66.847 


0.99943 


SDSSJ 005008.21 - 000359.0 


4 





24.004 


0.07446 


SDSSJ005245.il - 005337.2 


6 


5 


2127.066 


0.99999 



Notes. Name of the system, total number of spectra used for 
identification, number of spectra from own spectroscopic obser- 
vations, time span between first and last observation and prob- 
ability of measuring large radial velocity variations. 



4.2.1. Orbital periods from photometric observations 

Photometric variability in PCEBs can be caused by four 
different mechanisms: ellipsoidal modulation due to the dis- 
torted shape of the secondary star, reflection effect on the 
surface of the secondary star due to irradiation from the 
hot primary, eclipses, and rotational modulation due to 
starspots (BY Dra phenomenon). The first three types of 
variability are, by definition, modulated on the orbital pe- 
riod. For most PCEBs, also variability due to star spots 
will be modulated on the orbital period, as the companion 
will by tidally synchronised. However, for wide binaries, 
star spots on the companion star may mimic variability of 
a short-period PCEB. Four out of the seven WDMS bina- 
ries for which we obtained time-series photometry (Table[l} 
show ellipsoidal modulation, and we determined their or- 
bital periods by fitting sine curves to the observations. 
Periodograms and phase folded light curves are presented 
and discussed in more detail in Appendix [A"l 



The distribution of P(x 2 ) is shown in Fig. [5J again for 
the SDSS data (solid line), our follow-up data (dashed line), 
and the combined data (gray). The combined probability 
distribution presents a pronounced peak at P(x 2 ) = 1 
corresponding to systems exhibiting a statistically signif- 
icant radial velocity variation, i. e. PCEBs, and a less pro- 
nounced one close to P(x 2 ) = 0, corresponding mainly to 
very wide systems, and a small contribution of PCEBs with 
very low inclinations resulting in projected radial velocities 
that are too low to be picked up by the spectroscopic data. 
In total we found 191 PCEBs among the 1246 SDSS WDMS 
binaries for which at least two radial velocity measurements 
are available. These numbers can be used to set a strict 
lower limit on the fraction of PCEBs, see Sect. l6.1l for a more 
detailed analysis. For clarity, as for Fig. [TJ we mention that 
the combined distribution of P(x 2 ) is not the direct sum 
of the individual sets and illustrate this with an example: 
SDSS 0848+0058 has three SDSS subspectra taken within 
1.4 hours, and we took two spectra separated by five nights. 
We obtained P{\ 2 ) = 0.4573 and P{ X 2 ) = 0.4014 respec- 
tively. Combining the data we have 5 spectra, separated by 
7 years and obtain P(x 2 ) = 0.7140. 



4.2.2. Orbital periods from spectroscopic observations 

We obtained further spectroscopy for 68 of the 191 PCEBs 
identified by our radial velocity survey (Table[lJ. A period 
search was performed b y computing periodograms adopting 
ISchwarzenberg-Czernvt s (1996) ORT/TSA method within 
MIDAS, which fits periodic orthogonal polynomials to the 
observations, combined with an analysis of variance statis- 
tic. In a second step, we fitted sine curves of the form 



7scc + 



Pnvh 



(1) 



to the radial velocity data, with 7 SCC the systemic velocity of 
the secondary star, to the zero point defined by the inferior 
conjunction of the secondary star, and K sec the radial veloc- 
ity amplitude of the secondary star. Initial values for P or b 
were selected from the highest peaks in the periodograms. 
We adopt as orbital period the values of the best-fit solu- 
tion, which was always coinciding with the highest peak in 
the periodogram. We could measure orbital periods for 58 
of the 68 PCEBs that were spectroscopically followed up. 
For the remaining 10 systems the radial velocity data was 
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insufficient to unambiguously identify the true orbital pe- 
riods. The radial velocity variations measured for these 10 
systems range from 50 kms -1 to 490 kms -1 so a priori 
there is no reason to think that their orbital period dis- 
tribution differs significantly from the one drawn from the 
rest of the population of PCEBs. 

The 58 orbital periods range from 1.97 to 104.56 hours 
(0.08 to 4.35 days) (Table©. The spectra of these PCEBs, 
and their phase-folded radial velocity curves are shown in 
Fig. [3j [5] and Fig.|4l[6l respectively, and their stellar and bi- 
nary parameters given in Table [SJ which also lists K sec and 
7 SCC . The orbital period distribution of the SDSS PCEBs is 
discussed in more detail in Sect. [5] 



4.2.3. Orbital inclinations and candidate eclipsing systems 

We estimated the orbital inclinations and binary sepa- 
rations of those PCEBs for which masses of the white 
dwarf and the companion star (Sect. [2]), Ki and P or b 
were available using Kepler's third law, and propagated 
the errors in all quantities accordingly. We further calcu- 
lated if wd, as well as the Roch e-lobe size of the com- 
panion star, i?i b using lEggletonl 's (1983) approximation. 
Fi nally, using the empirical radiu s-spectral type relation 
of iRebassa-Mansereas et all <|2007t) we estimate the Roche- 
lobe filling factors of the secondary stars. All this ancillary 
information is given in Table[5j In case the effective temper- 
ature of the primary is smaller than 12 000 K, the derived 
spect r oscopic mass is inher ently uncertain ()Koester et al.l 
2009; iTremblav et al.l I2010D and the derived parameters 
should be used with caution. 

For six systems, the estimated orbital inclinations are 
close to 90 degrees, making them good candidates for be- 
ing eclipsing binaries: SDSS 0238-0005, SDSS 1348+1834, 
SDSS 1434+5335, SDSS 1506-0120, SDSS 2132+0031 and 
SDSS 2318-0935. Eclipsing systems offer the best oppor- 
tunity to derive accurate fundamental stellar parameters, 
i.e. masses and radii, and we have begun to exploit the 
potential offered by the eclipsing PCEBs (ISteinfadt et all 
120081: iNebot Gomez-Moran et aljfcdbH iPvrzas et al.ll2009D . 
SDSS 1348+1834, aka CSS 21357, was recent ly identified 
in th e Catalina Real Time Transient Survey (jDrake et al.l 
120091) as an eclipsing PCEB (jDrake et al.ll2010D . confirming 
the predictive power of our estimates. 



4.2.4. Observed orbital period distribution 

We have measured here the orbital periods of 58 PCEBs. 
Together with the 15 period s already published by our 
group (IPvrzas et al. 120111: pRebassa-Mansergas et al.l 



0.1 



Portrays] 



2008: iSchreiber et al.1 120081: IPvrzas et all 12 009: 
Nebot Gomez-Moran et al.ll2009l) . two more bv lDrake et al.l 
(120101) . and four additional PCEBs that are subject to 
forthcoming individual studies (Pyrzas et al in prep.; 
Parsons et al in prep.), the total number of PCEBs from 
SDSS with orbital period measurements stands at 79, 
which is more than double the corresponding number of 
pre-SDSS PCEBs. The orbital period distribution of the 
SDSS PCEBs is shown in Fig. [7] on a logarithmic scale, in 
light and dark gray those systems that were identified as 
PCEBs based on spectroscopy (separating systems with 
PCEB identification spectra taken in the same night, light 
gray) and photometrically identified PCEBs are shown in 
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Fig. 7. Orbital period distribution of the 79 SDSS PCEBs 
characterised to date. Systems identified as PCEBs based 
on spectroscopic observations are shown in light gray if the 
spectra were taken during the same night and dark gray 
otherwise. In medium gray we show those systems that were 
identified as PCEBs based on photometry. The 58 systems 
added by this paper to the orbital period distribution are 
indicated with vertical lines. 



medium gray. The 58 periods determined in this paper are 
indicated by vertical lines. Taking this observed period 
distribution at face value, PCEBs are found to have 
orbital periods ranging from 1.97 h to 4.35 d, they follow 
approximately a normal distribution in log(P or b) with a 
peak centred on ~ 8.1 h and a width of 0.41 in log(P rb[d]) 
as determined from a Gaussian fit. 



5. Analysis 

We have performed a radial velocity survey of 1246 WDMS 
binaries. The majority of the systems have three radial 
velocity measurements, spread over hours to hundreds of 
days (see Fig. [I}. We have identified 191 of these systems as 
PCEBs based on the detection of significant radial velocity 
variations. While we consider the remaining 1055 systems 
to be candidate wide binaries, it is clear that there will 
still be a number of PCEBs among them which our radial 
velocity survey failed to identify because of observational 
selection effects. 

Our ultimate aim is to establish the intrinsic period dis- 
tribution of the SDSS PCEBs. We will first investigate in 
the following sections what information might be gleaned 
from the distribution of probabilities of our initial radial 
velocity survey (Fig. [3]). Subsequently, we will analyse the 
observational biases within the orbital period distribution 
of the SDSS PCEBs (Fig.[7J|, and derive a bias-corrected 
period distribution, that will then be discussed in Sect. El in 
the context of compact binary evolution. 

5.1. Does our radial velocity survey constrain the orbital 
period distribution of WDMS binaries? 

The detailed shape of the probability distribution used to 
identify PCEBs among SDSS WDMS binaries will depend 
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Fig. 8. Assumed orbital period distributions (left panels) 
and associated simulated distribution of the probability of 
measuring large radial velocity variations (gray) with mea- 
sured distribution in black for comparison (right panels). 




l°g(P ort ,[days]) 

Fig. 9. Detection probability for three cases corresponding 
to 2 spectra taken at the VLT (green), 6 spectra taken at 
the WHT, Magellan and SDSS (blue) and 16 spectra from 
the SDSS (red). 



on the technical details of our radial velocity survey (spec- 
tral resolution, temporal sampling, number of spectra) and 
the intrinsic period distribution of the WDMS binaries. To 
test how sensitive the observed probability distribution is 
to differences in the underlying orbital period distribution, 
we performed extensive Monte Carlo simulations. As in- 
put, we use the measured stellar masses (Table©, the times 
of observation, and the radial velocities and associated er- 
rors (Table© for all 1246 WDMS binaries within our ra- 
dial velocity survey. We then adopted four different orbital 
period distributions: (1) short-period binaries only, with 
2h < P or b < 100 d, and uniformly distributed in log(P or b); 
(2) wide binaries only, with 3.5 < log(P or b[d]) < 6.5, uni- 
formly distributed in log(P or b); (3) a bi-modal distribution 
with a short-period binary population peaking at 7.5 h, nor- 
mally distributed in log(P or b) and with a minimum period 
of ~ 1.5 h, plus a second population of wide binaries with 
3.5 < log(P or b[d]) < 7.5, uniformly distributed in log(P or b), 
and finally (4) a uniform period distribution in log(P or b) 
with log(P orb [d]) < 8.5. Obviously (1), (2) and (4) are very 
unlikely to be physically relevant, but were adopted to run 
the test with several radically different period distributions. 
We then used 10 4 random inclinations per object, uniformly 
distributed in cosi, and calculated the probability distribu- 
tions of the two adopted period distributions. We assumed 
circular orbits for simplicity (and note that the orbits of 
PCEBs will be, in most cases, circularised). 

The results of this simulation are illustrated in Fig. [8j 
showing the adopted period distributions on the left hand 
side, and the resulting probability distributions on the right 
hand side, with the observed probability distribution shown 
for comparison. The top two panels show that a peak at 
P(x 2 ) = 1 requires close binaries, whereas a broad accu- 
mulation towards P(x 2 ) — is dominated by wide bina- 
ries. In other words, the observed probability distribution 
(Fig.© can only be explained if the underlying orbital pe- 
riod distribution of the observed WDMS binaries contains 
both wide binaries (P or b~10 3 d) and PCEBs with shorter 
periods (P orb < lOOd). 



However, the bottom two panels show that the two pe- 
riod distributions which contain both short-period PCEBs 
and wide binaries, but radically differ in shape, both lead 
to probability distributions that are very similar to the ob- 
served one. We conclude that our radial velocity survey car- 
ried out to identify PCEBs alone cannot constrain the de- 
tailed shape of the underlying WDMS binary orbital period 
distribution, and that it is hence necessary to accurately 
measure the orbital periods for a large and homogeneous 
sample of PCEBs. 



5.2. Observational biases and selection effects 



We have presented in Sect. l4.2.4l the orbital period distri- 
bution of 79 PCEBs that were spectroscopically identified 
by SDSS (Fig.©. The resulting orbital period distribution 
peaks at ~ 8 hrs and very few systems with orbital periods 
longer than a day have been found. 

This is not only the largest but also the most ho- 
mogeneous sample of PCEBs available so far, which 
has a large potential to provide crucial constraints on 
the theories of close compa c t bin a ry formation and evo- 



lution (e.g. ISchreiber et all l2010t IZorotovic etall 120101: 



iRebassa-Mansergas et al.l 201 1|) . It is, however, of funda- 
mental importance to understand any observational selec- 
tion effects that may affect the properties of the SDSS 
PCEB sample. In the context of the orbital period distri- 
bution, there are two particular points that need to be in- 
vestigated. Firstly, our radial velocity PCEB identification 
method is obviously more sensitive to short orbital period 
and high inclination systems and secondly, we only mea- 
sured the orbital periods of a sub-sample of the identified 
PCEBs which is not necessarily representative for the entire 
sample. 
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5.2.1. PCEB identification 

While we have compiled the largest and most homogeneous 
sample of PCEBs to date, not all the systems in this sample 
have been identified in the same way. 

Ten systems were initially classified as PCEBs because 
of photometric variability and four because of radial ve- 
locity variations seen over the course of a single nighfl 
Both methods favour the identification of short orbital pe- 
riod PCEBs. This is illustrated in Fig. where all those 
14 systems have periods Ss 9h (0.4 d), and we discard them 
in the following analysis. 

The remaining 65 PCEBs were identified by radial 
velocity variations from spectra that have been obtained on 
different nights, and we will concentrate our analysis of the 
observational biases on this "core" sample. The strength of 
the observational bias introduced in a radial velocity search 
for binaries depends on the intrinsic orbital period and the 
sampling frequency. 

We have performed Monte -Carlo simulatio ns sim- 
ilar to those presente d in ISchreiber et al. 
iRebassa-Mansereas et al.l ([20081 ): 



200S 

Schreiber et~aT1 (1201C 



to evaluate if the observed orbital period distribution 
(Sect. 14.2.21) is significantly biased by the design of our ra- 
dial velocity survey, and we calculated the detection proba- 
bility as a function of orbital period for all 65 PCEBs. Given 
the times of observation and the error of the radial velocity 
measurements, we randomly select 10 4 phases and inclina- 
tions for a given orbital period, and determine the corres- 
ponding x 2 probability. The PCEB detection probability 
for each system as a function of the orbital period is then 
simply given by the number of cases with P(x 2 ) > 0.9973 
(> 3ct) divided by 10 4 . The resulting detection probabili- 
ties depend on the number of spectra taken, the temporal 
sampling of the measurements, and the accuracy of the ra- 
dial velocity measurement (i.e. basically the spectral res- 
olution of the telescope/instrument combination that has 
been used). 

Figure [5] shows the obtained detection probabilities for 
three rather extreme examples: SDSS 0807+0724 (green, 
P orb = 11.45 h, identified as a PCEB from two VLT/FORS 
spectra), SDSS 0152-0058 (blue, P orb = 2.15 h, 3xWHT, 
2 x Magellan, lxSDSS), and SDSS 0225+0055 (red, P orb = 
21.86 h, 16 SDSS subspectra). We have chosen these three 
cases to illustrate that our PCEB identification method was 
sensitive to PCEBs with orbital periods in the range of 1 — 
10 days (detection probabilities P(x 2 )~0.35) as long as the 
time baseline exceeded one night. In addition, Fig. |H] illus- 
trates the following: (1) the lower the resolution, the lower 
is the detection probability, even in the most favourable 
case, i.e. a large number of spectra. (2) Variations in the 
detection probability with orbital period are dominated by 
phase coverage and are therefore strongest when only few 
spectra are available. Orbital periods near one day and their 
first subharmonics (12 h, 8h) have substantially depressed 
detection probabilities, as this increases the chances to sam- 
ple the same orbital phase twice on two subsequent nights. 
(3) Superimposed on the fine structure related to the details 
of the sampling is a decreasing envelope, which is caused by 
a larger range of (low) inclinations leading to radial velocity 



6 SDSS 0106-0014, 0110+1326, 0303+0054, 0908+0604, 1212- 
0123, 1423+2409, 1548+4057, 1648+2811, 1724+5620 and 
2112+1014 

7 SDSS 0307+3848, 1300+1908, 1348+1834 and 1625+6400 
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Fig. 10. Orbital period distribution of the 65 SDSS PCEBs 
discovered in an homogeneous way through radial veloc- 
ity variations measured in spectra spread over different 
nights (gray) and the bias-corrected distribution (white). 
The Gaussian fit to the bias-corrected distribution is plot- 
ted with a dashed line. In the upper panel we show the 
mean probability of detection as a function of the orbital 
period (see text) and the corresponding mean value at each 
bin centre (red dots). 



variations that are too low to be resolved. The mean detec- 
tion probability for the whole sample, calculated from the 
65 individual detection probabilities, is higher than 35% up 
to 10 days orbital periods (Fig.[TUl upper panel), demon- 
strating that our radial velocity survey is well capable of 
identifying PCEBs with long orbital periods. We will use 
this mean detection probability to bias-correct the observed 
orbital period distribution in Sect. 15.2.31 

5.2.2. PCEBs selected for orbital period measurements 

Given that we have measured orbital periods only for a frac- 
tion of the PCEBs that we identified, the next question to 
address is whether the subset of systems is representative 
of the whole sample. For instance, our follow-up strategy 
might have been biased towards the "easier" targets, i.e. 
those which exhibited large variations in our radial veloc- 
ity survey. For this purpose, we compared the cumulative 
radial velocity variations of the 65 PCEBs identified in our 
radial velocity survey (see Sect 14. ip . and of those PCEBs 
for which we obtained no additional follow-up observations 
(Fig. [11]). A Kolmogorov-Smirnov test between both sets 
gives a probability of 0.97 of both samples being drawn 
from the same parent population. We can conclude that 
both distributions are very similar, and that we have mea- 
sured the orbital period of a representative sub-sample of 
the SDSS PCEB population. 

5.2.3. Bias corrected orbital period distribution 



We have demonstrated in Sect. 15. 2.21 that the 65 PCEBs 
that were identified from multiple-night radial velocity 
snapshots and have measured orbital periods are represen- 
tative of the entire SDSS PCEB sample. Here, we apply 
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Fig. 11. Cumulative distribution function (CDF) of the ra- 
dial velocity variations used for the identification of PCEBs, 
systems with measured orbital period are shown in red and 
systems that were not further followed-up are shown in 
black. The gray area is the corresponding CDF taking the 
la radial velocity errors into account). 



Two observational selection effects will affect this em- 
pirical estimate of the PCEB fraction. On the one hand, 
our analysis in Sect . 15 . 2 .T1 showed that a number of PCEBs 
will escape identification because of their low orbital incli- 
nation. From the mean detection probability and the bias 
corrected orbital period distribution we learn that ~ 16% 
of the PCEBs are not detected due to a combination of in- 
clination and phase coverage (see Fig. [TU)) . Assuming that 
the 164 detected PCEBs (observed in different nights) fol- 
low the same orbital period distribution as the observed one 
(drawn from 65 of them) the number of expected PCEBs 
is 195, implying a fraction of PCEB to WDMS binaries 
of 27 ± 2% (195/718). On the other hand, WDMS bina- 
ries that are sufficiently wide will appear as blended ob- 
jects, or even fully resolved objects, and not have been fol- 
lowed up spectroscopically, as they have unreliable colours 
or colours of single white dwarfs and M-dwarfs, respectively. 
Hence, our estimate of 27 ± 2% PCEBs among the spa- 
tially unresolved SDSS WDMS binaries can be taken as 
an upper li mit on the intrinsic fra ction of PCEBs among 
all WDMS. ISchreiber et al.l (|2010ft estimated the fraction 
of SDSS resolved WDMS to be 12-23%. Taking this into 
account we finally estimate the intrinsic fraction of PCEBs 
to be - 21 - 24%. 



the average detection probability fSect. l5T2~Tj) to correct 
their observed period distribution (dark gray histogram in 
Fig. ITU1) . The upper panel of Fig.fTU] shows the detection 
probability averaged across all 65 systems, and its binned 
values as red dots. The bias-corrected period distribution 
is shown in Fig. [TU] (white histogram). It is evident that the 
observational bias of our radial velocity survey does not 
dramatically affect the conclusions drawn in Sect. 14 .2.11 i.e. 
the period distribution of the SDSS PCEB sample follows 
broadly a normal distribution centred at ~ 10.3 h (0.43d), 
a width of 0.44 in log(P rbM]), a short-period cut-off at 
~ 2 h and a long-period cut-off at ~ 5 d, where the position 
of the peak and the width have been determined from a 
Gaussian fit (dashed line). The bias-corrected period dis- 
tribution is discussed in the context of published binary 
population models in Sect. 16.41 

6. Discussion 

6.1. The PCEB fraction among WDMS binaries 

Analysing the radial velocity variations of 1246 SDSS 
WDMS binaries we identified 191 PCEBs (Sect. ED}, which 
provides a lower limit of ~ 15±1% to the fraction of PCEBs 
among the SDSS WDMS binaries. However, for a substan- 
tial number of these 1246 systems the only available source 
of radial velocities are the SDSS subspectra, which, as they 
are usually taken during a single night (see FigfT]), are sen- 
sitive only to very short orbital periods. 

If we only consider the sub-sample of 718 WDMS bi- 
naries that were observed on different nights, we find 164 
PCEBs, i.e. a PCEB frac tion of ~ 23 ± 2% This number is 
similar to that found by ISchreiber et all (|2008l ). based on 
our VLT/FORS radial velocity survey of 26 SDSS WDMS 
binaries. This empirical determination of the PCEB frac- 
tion is broadly consistent with t he results of the binary p op- 
ulation synthesis carried out by IWillems fc Kolbl (2004) . 



iFarihi et"a l. (2006, 201(|) pursued a complementary ap- 
proach to investigate the make-up of the known WDMS 
binary population, obtaining high spatial resolution imag- 
ing with HST/ ACS. They observed a total of 72 WDMS 
binaries unresolved from the ground, and found a bimodal 
distribution in projected binary separations: all spatially 
resolved systems had binary separations J> 10 AU, clearly 
above the detection threshold of 1 — 2 AU. Accounting for 
the freque ncy of WDMS binaries th at can be resolved from 
the ground lFarihi et al.l (|2005l . l2010D conclude that the frac- 
tion of short-period binaries among all WDMS binaries is 
~ 25%, which is in excellent agreement with our estimate 
presented here. 19 of the WDMS binaries imaged with HST 
were also spectroscopically re-discovered by SDSS, and are 
hence in our WDMS binary sample (Table[6]). 

IFarihi et all ()2010D predict that all the WDMS bina- 
ries unresolved by HST (Table 5 in their paper) should be 
short-period PCEBs, which is confirmed for all five systems 
for which we have obtained follow-up observation, span- 
ning more than one night. For an additional two of Farihi's 
unresolved systems (SDSS1341+6026, 1500+1659) we have 
only the radial velocity information of the SDSS subspectra 
spanning less than one hour, and we classified both systems 
as wide binary candidates. This underlines what we said 
above, i.e. that the single-night SDSS subspectra allow us 
to identify only the shortest-period WDMS binaries (which 
we fully take into account by our Monte Carlo analysis of 
the selection biases). For complet eness, we no t e that one 
additional unresolved binary from IFarihi et al.l (|2010D was 
recently found to be a n eclipsing W DMS binary with an 
orbital period of 2.93h lDrake et all (|2010D . 

Vice versa, based on our radial velocity information 
for nine of the WDMS binar ies spatially resolved by HST 
(Table 4 of IFarihi et aTll2010h . we classified all of them as 
wide binary candidates (labeled as "wide" in the bottom 
portion of Table [5]) . 
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Table 6. SPS S- WDMS binaries in common with 
iFarihi et al.l (|2010h . 



WD-Name SDSS-Name Type # spec At [d] 



Not resolved by HST 



WD0303 


- 007 


SDSS J030607.18 


- 003114.4 


pceb 


4 


2991 


WD0908 


+ 226 


SDSS J091143.09 


+ 222748.8 


pceb 


4 


5 


WD1051 


+ 516 


SDSS J105421.97 


+ 512254.2 


pceb 


3 


1089 


WD1140 


+ 004 


SDSS J114312.57 


+ 000926.5 1 


pceb 


3 


2551 


WD1339 


+ 606 


SDSS J134100.03 


+ 602610.4 


wide 


3 


0.03 


WD 1433 


+ 538 


SDSS J143443.24 


+ 533521.2 


pceb 


3 


1157 


WD1458 


+ 171 


SDSS J150019.33 


+ 165914.4 


wide 


2 


0.01 


WD1522 


+ 508 


SDSS J152425.21 


+ 504009.7 2 


pceb 


3 


1007 






Resolved by HST 








WD0257 


-005 


SDSS J030024.57 


- 002342.0 


wide 


2 


10 


WD1049 


+ 103 


SDSS J105227.72 


+ 100337.6 


wide 


3 


0.03 


WD1106 


+ 316 


SDSS J110843.03 


+ 312356.1 


wide 


3 


0.03 


WD1108 


+ 325 


SDSS J111045.90 


+ 321447.2 










WD1133 


+ 489 


SDSS J113609.59 


+ 484318.9 










WD1218 


+ 497 


SDSS J122105.34 


+ 492720.5 


wide 


3 


0.03 


WD1236 


- 004 


SDSS J123836.34 


- 004042.2 


wide 


i 


3197 


WD1419 


+ 576 


SDSS J142105.31 


+ 572457.1 


wide 


i 


0.05 


WD 1443 


+ 336 
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Notes. "Type" is our classification based on the radial velocity 
information. We list the number of spectra used for the classi- 
fication and the time span between them. Two of the systems 
unresolved by HST only have SDSS subspectra taken within 
< 1 h. 1 P oth = 9.274 h. 2 P orb = 14.15 h. 

6.2. PCEBs and detached CVs in the orbital period gap 

The observed orbital period distribution of CVs presents 
a sharp decrease of systems between 2 -3 h (Wa rnerlll976t 
IWhvte fc Eggletonlll980l: lKniggd[2006h . referred to as the 
orbital period gap of CVs, which is interpreted as a key to 
understanding their long-term evolution. The mass ratio of 
CVs implies that they are undergoing stable mass trans- 
fer, and it is angular momentum loss that drives them to- 
wards shorter orbital periods, keeping the companion stars 
in Roche-lobe contact. Above the period gap, the secondary 
stars have a radiative core, and are subject to magnetic 
wind braking that acts as a strong agent of angular momen- 
tum loss. The corresponding mass loss rates are sufficient to 
drive the companion stars of these long-period CVs out of 
thermal equilibrium. At ~ 3 h, the companion stars become 
fully convective, and magnetic wind braking is thought to 
stop. From there on, the evolution towards shorter periods 
is driven by gravitational wave radiation, resulting in much 
lower mass loss rates. Consequently, the secondary stars 
have time to relax to their thermal equilibrium radii, shrink 
within their Roche lobes, and turn into detached WDMS 
binaries with no mass transfer. Gravitational wave radia- 
tion keeps driving these detached systems towards shorter 
periods. At ~ 2h the secondary stars re-establish Roche- 
lobe contact, mass transfer stars again, and the systems 
re-appear as CVs. This "disrupted magnetic braking" sce- 
nario is currently the s tandard paradigm of CV evolution 
(jRappaport et all 1 19831) . 

Observationally, there is very little difference between 
genuine pre-CVs that will start mass transfer within, or 
shortly below, the CV period gap, and those detached CVs 
(dCVs) that are crossing the gap. The common hallmark of 
both populations is that they are WDMS binaries with sec- 
ondary spectral types ^M6 to ^M4. A very simple, but so 
far untested, prediction of the disrupted magnetic braking 



scenario is that there should be an excess of systems with 
periods of 2-3 h in the period distribution of WDMS bina- 
ries, due to dCVs that are crossing the orbital period gap. 
iDavis et~a l. (2008) carried out a number of binary popu- 
lation studies to estimate the strength of this effect, and 
predict that within the CV period gap the number of dCVs 
should exceed that of genuine pre-CVs with secondary stars 
in the mass range 0.17 — 0.36M Q by a factor ~ 4 — 13, which 
implies an increase of systems (PCEBs plus dCVs) of 1.5- 
2 in the gap in comparison with PCEBs at longer periods 
(see their Figs. 4 &10). 

The large sample of SDSS PCEBs provides a first op- 
portunity to test this prediction. Observationally, the sec- 
ondary stars in the SDSS WDMS binaries are characterised 
in terms of their spectral type. CVs just above and below 
the period gap have ~M3-4, and ^M5 secondary stars, re- 
spectively. We show in Fig.Q^the period distributions of all 
SDSS WDMS binaries with M3-4 companions (top panel), 
and a somewhat broader range of M2-5 (bottom panel) 
which accounts, somewhat conservatively, for the uncer- 
tainty in the spectral types of the companion stars. The 
expected increase towards shorter periods by a factor of 
1.5-2 is not present in our data. In contrast the number of 
dCVs plus PCEBs is slightly decreasing (M3-4) or constant 
(M2-5), although the number of systems in this sub-sample 
is still rather small. Comparing the number of dCVs plus 
PCEBs in the gap with the number of PCEBs between 3-4 
hours, we obtain a fraction of dCVs plus PCEBs (in the 
gap) to PCEBs (above the gap) of 0.25 ± 0.30 (M3-4) or 
0.43 ± 0.30 (M2-5), which excludes the predicted increase 
of 1.5-2 on a 3ct (M3-4) or a 2cr (M2-5) level. 

Three SDSS PCEBs that match the expected characte- 
ristics of dCVs are SDSS 0052-0053, SDSS 1611+4640 and, 
SDSS 2243+3122, with orbital periods 2.74, 1.98 and 2.87h 
respectively. The Roche-lobe filling factors of SDSS 0052- 
0053 and SDSS 1611+4640 are 1.0 ± 0.6 and 1.52 ± 0.50 
respectively. For SDSS 2243+3122 we have no estimate of 
the (DC) white dwarf mass, and hence can not determine 
the filling factor. 

6.3. Comparison with non SDSS PCEBs 

In this section we compare the orbital period distribution 
of PCEBs discovered within the SDSS, with that of PCEBs 
that were discovered by other surveys. We compiled a list of 
48 non-SDSS PCE Bs containing a white dwarf plus a main 
sequence star fromlRitter fc Kolbl (|2003l version v7.14) and 
iBrown et aD (|201ltl . The orbital period distribution of non- 
SDSS PCEBs (Figure IT3l gray histogram) is overall similar 
to that of the SDSS PCEBs (white hatched histogram) , pre- 
senting both a peak at around ~ 5 h and decreasing number 
of systems towards long orbital periods. A Kolmogorov- 
Smirnov test between the cumulative period distributions 
(bottom panel in Figure [T2]) gives a probability of 35%, i.e. 
consistent with the two samples being drawn from the same 
parent population. 

We not e that this resu l t is in agreement with the con- 
clusions of Miszal ski et al.l (|2009f) , who first noticed the re- 
semblance between the orbital period distribution of binary 
central stars in planetar y nebulae (PNe) and that one of 
PCEBs as presented bv iRebassa-Mansergas et~a l. (2008). 
This was later confirmed by iHillwid (|201l[ ) , who selected 
27 WDMS binaries that are the central stars of plane- 
tary nebulae (PNe) and the 35 SDSS WDMS binaries (as 
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Fig. 12. Orbital period distribution of PCEBs with sec- 
ondary spectral types in the range M3-M4 (upper panel) 
and M2-M5 (bottom panel). Magnetic braking disruption 
predicts a peak at around 2-3 hours corresponding to an 
increase in systems due to detached CVs, we highlighted 
the region in gray. 




P„,b Ways] 

Fig. 13. Orbital period distribution of the 48 non-SDSS 
PCEBs that contain a white dwarf plus a main sequence 
star (gray histogram) and the 65 systems SDSS PCEBs 
discovered through radial velocity variations on spectra 
taken on different nights (dashed histogram) in days. 
Bottom panel: cumulative distribution function of non 
SDSS PCEBs (gray line) and SDSS PCEBs (black dashed 
line). 



published bv lZorotovic et all (|2010t) V Surveys dedicated to 
measure the orbital period of the binar y central sta rs of 
PNe have been mostly photometric (see iBondl l2000l and 
referenc es therein), implying a bias towards short orbital 
periods. iDe Marco et al.l ((2008) studied the possible biases 
introduced in the photometric sample, and concluded that 
the decrease of the number of PNe towards long periods is 
unlikely to be due to observational biases. 



6.4. Comparison with predictions from BPS models 

In this section we discuss the SDSS WDMS binary sam- 
ple in the context of published binary population synthe- 
sis models. However, it is necesary to keep in mind that 
the SDSS PCEB sample is subject to selection effects, i.e. 
systems with companion spectral types earlier than ^M0 
are severely underrepresented, and a comparison of the pre- 
dicted orbital period distribution of the different BPS stud- 
ies for the corresponding mass domain is in most cases not 
possible. 

Common to all BPS models is the prediction that the 
population of WDMS binaries should be made up of wide 
binaries, with orbital periods longer than one year, in which 
the components never interacted, and close binaries which 
formed from common envelope evolution. The two most 
important parameters determining the outcome of the CE 
phase are the initial mass ratio distribution of the progeni- 
tor binaries, and the efficiency of the ejection of the common 
envelope, acE, which will affect the position of the peak in 
the orbital period distribution of the PCEBs, and the slope 
towards longer orbital periods. 

On e of the f i rst B PS fo r WDMS binaries was c arried 
out bv Ide Kooll ([1992) and Ide Kool fe Ritterl (|l993l) . pre- 
dicting that most PCEBs should have periods shorter than 
~ one day. An extended tail towards longer orbital peri- 
ods up to 1000 days is expected for «ce = 1, while no 
tail is observed for «ce = 0.3. The authors also simulated 
the expected properties of a volume limited blue sample 
of present-day WDMS binaries, however, comparison with 
the very small number of PCEBs at that time did not 
lead to any significant observational constraints. Our or- 
bital period distribution (Fig. ITUl) is very similar to their 
predicted one for acE = 0.3 and the mass of the secondary 
sta r picked randomly fro m the IMF, already pointed out 
by iMiszalski et al.l (|2009l ) in the context of PN, or with 
dN tx dq (see their Fig. 12b). However , this comparison 
should be regarded with some caution, as ide Kool fe Ritterl 
(|l993l ) simulated their "observational" sampl e for the se- 
lectio n criteria of the Palomar- Green Survey (|Green et al.l 
1986) (U < 16 and U - B < -0.5), whereas the SDSS 
WDMS sample has a much larger limiting magnitude i ~ 
19.1 and is drawn from a much larger ugriz colour volume. 

IWillems fc Kolbl (|2004l ) carried out an extensive BPS 
study of WDMS binaries, using acE = 1, and for different 
initial mass ratio distributions. Their BPS predicts that the 
orbital period distribution of PCEBs peaks around 1 day, 
with a substantial tail towards longer orbital periods of up 
to 100 days, which is in contrast to the p eriod distribution 
of the SDSS PCEBs (Fig. [p. However. IWillems fc Kolbl 
(|2004l ) did not convolve their model WDMS binary popu- 
lation with detailed observational selection effects, but just 
discussed their results in terms of luminosity ratio distri- 
butions, iwd/^MSj which prevents any more quantitative 
comparison with our SDSS WDMS binary sample. 

iPolitano fc Weilerl (|2007l ) investigated the effects on the 
population of PCEBs for ctcE depending on the mass of the 
secondary star. They predict an orbital period distribution 
peaking at around 3 days, which is incompatible with the 
observed peri od distribution of the SD SS WDMS binaries. 
However, also IPolitano fc Weilerl (|2007l) did not simulate a 
realistic "observational" WDMS population, again prevent- 
ing a quantitative comparison with our sample. 
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The most recent BPS o f WDMS binaries was per- 
formed bv I Davis et al.l (|2010f ). who, in contrast with previ- 
ous works, calculate the binding energy parameter A from 
stellar evolution models incorporating the internal energy 
of the envelope. While this is a significant step forward com- 
pared Tto_p2evious PBS that worked with a constant value 
of A, [Davis et al.l (12010T) o btain ed the values of A from the 
tables of iDewi fc Taurisl ((2000) , which cover only a very 
small part of the re levant parameter sp ace. The orbital pe- 
riod distribution of iDavis et al.l (|2010l) predicts a long tail 
of systems with periods > 1 d, inconsistent with the period 
distribution of the PCEBs known at the time of publication 
(their Fig. 10), even when taking into account th e detection 
probabilities of lRebassa-Mansergas et all ()2008l ) . Our more 
detailed observed (FigJT]) and bias-corrected (Fig.[TU]) or- 
bital period distributions are in better agreement with their 
predicted orbital period distribution limited to PCEBs with 
companion masses < 0.5 Mq, i.e. spectral types later than 
^M0 (see their Fig. 10 middle panel). 

In summary, most published BPS models predict a sub- 
stantial number of PCEBs with orbital periods > 1 d, which 
is currently incompatible with the orbital period distribu- 
tion of the SDSS PCEBs. The observed dearth of long- 
perio d PCEBs could be re conciled with a low value of 
acE- IZorotovic et al.l (|2010f) used a sub-set of 35 of the 
PCEBs presented here, along with some other 25 non-SDSS 
PCEBs, to constrain the value of «ce by reconstructing 
their evolution. They found that there seems to be no de- 
pendence of acE on the mass of the secondary star, and 
that if the value is universal only a small value of acE 
can explain the short orbital periods observed among the 
known PCEBs. However, it is also clear that even the SDSS 
PCEB sample is subject to selection effect, i.e. systems with 
companion spectral types earlier than ^M0, i.e. masses 
> O.5M are severely underrepresented, which are ex- 
pected to have longer orbital periods than the currently 
dominant population of PCEBs with low-mass (mid M- 
type) companions. 

7. Conclusions 

Nearly 2000 WDMS binaries have been spectroscopically 
identified within the SDSS survey. A detailed characterisa- 
tion and analysis of this large and homogeneously selected 
sample promises to improve substantially our understand- 
ing of close binary evolution, in particular of the common 
envelope phase. Using follow-up observations of 385 WDMS 
binaries, combined with the SDSS subspectra for an addi- 
tional 864 systems, we identified 191 PCEBs. We determine 
orbital periods for 58 of these systems, increasing the total 
number of SDSS PCEBs with known orbital periods to 79. 
We analyse the observational biases inherent to the SDSS 
PCEB sample, and find that 21 - 24% of all SDSS WDMS 
binaries with M dwarf companions must have undergone 
common envelope evolution. The bias-corrected orbital pe- 
riod distribution of the SDSS PCEBs shows a peak near 
~ 10.3 h, and drops fairly steeply to both shorter and longer 
periods, with no systems at P rt> < 1.9 h, and few systems 
at P > 1 d. Comparing the properties of the SDSS PCEB 
population to the results of published binary population, 
we conclude that the observed dearth of long-period sys- 
tems can be explained by low values of «ce- Observational 
selection effects against earlier-type companions can also 
play a role. Given a Roche lobe-filling giant or asymptotic 



giant (the progenitor of the white dwarf component) in a 
binary system, the more massive the companion star, the 
greater the orbital energy available to eject the common 
envelope, and therefore the longer the orbital period of the 
remnant PCEB, other things being equal. Unlocking the di- 
agnostic potential of the SDSS WDMS binary population 
now requires the computation of a state-of-the art binary 
population model that includes a realistic simulation of the 
observational selection effects inherent to this sample. 
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Appendix A: Optical light curves 

Photometric follow-up observations were carried out for a 
total of seven WDMS binaries where the available spec- 
troscopy suggested a short orbital period. Telescope, fil- 
ter and duration of the observations are listed in Table [T] 
Four out of the seven observed systems showed el- 
lipsoidal modulation: SDSS 0853+0720, SDSS 1611+4640, 
SDSS 2112+1014, and SDSS 2243+3122. We determined 
their orbital periods by computing periodograms and fit- 
ting sine curves to the phase folded light curves (see 
Fig. EU) . We estimated P orb = 3.6,1.9,2.2, and 2.8 
hours respectively. Spectroscopic follow-up observations 
were carried out for SDSS 0853+0720, SDSS 1611+4640, 
and SDSS 2243+3122 covering a full cycle, and confirm- 
ing the measured periods. We could connect the spec- 
troscopic and photometric runs for SDSS 1611+4640, and 
SDSS 2243+3122, without cycle count aliases, which al- 
lowed us to refine their orbital periods. 

The light curves of SDSS 0853+0720 and 
SDSS 2243+3122 show both two uneven maxima at 
phases 0.25 and 0.75. For SDSS 0853+0720 the first 
maximum is brighter than the second one, while for 
SDSS 2243+3122 is the other way round. In both systems 
equal minima at phases and 0.5 are observed and the 
maximum variation is almost 0.2 magnitudes in the r 
and the i bands respectively. Ellipsoidal modulation and 
spots in one of the hemispheres of the secondary star 
could explain the uneven maxima (O'Connell effect, see 
iLiu fc Yand (|2003h ). Photometry carried out on the 15th 
of August 2009 for SDSS 2243+3122 revealed a flare 
with ~ 25 minutes length in the decay, and rising by 
0.35 magnitudes in the i band (see Figure IA.2|) . This 
would strengthen the argument of magnetic activity being 
the cause of the observed disparity in the maxima (see 
Figure IA.1[) . 

Although no radial velocity curve is available for 
SDSS 2112+1014, the five available spectra show radial ve- 
locity variations of more than 300 kms -1 in less than 40 
minutes, indicating a rather short orbital period, consistent 
with the measured photometric orbital period of 2.2 hours. 
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Fig. A. 2. Differential photometry in the i band obtained 
for SDSS 2243+2122 the 15th August 2009. A flare of ~ 
25 minutes length with a relative enhancement of ~ 0.4 
magnitudes was observed. 



Appendix B: Notes on individual systems 

SDSS 0225+0054, the prim ary star was classified as a 
DC bv lKleinman et all (|2004l) . Although it is true that no 
absorption lines are visible in the primary star we con- 
sider that the signal to noise of the residual white dwarf 



after the subtraction of the best fitted secondary tem- 
plate is insufficien t for giving any other classifica tion than 
white dwarf (see Rebassa -Mansergas et alj 120071 for de- 
tails on the spectral analysis). SDSS 0225+0054 together 
with SDSS 0853+0720 are WDMS binaries consisting 
of cold DC white dwarf and an M4.5 and M3 compan- 
ion star, clear examp l es of the old PCEBs predicted by 
iSchreiber fe Gansickd (|2003f ) that would not have been 
found in any search for blue objects. 

SDSS 0307+3848 (SDSS J030716.44+384822.8) has a red 
star nearby - within 5"- (SDSS J030716. 16+384821.4) and 
is slightly brighter in the r band (16.70 versus 15.17). To 
be sure that there was no contamination in our spectra 
by this source we rotated the slit so as to take spectra si- 
multaneously and compared the radial velocities with those 
of this system. No variations in the radial velocity of this 
visual companion were observed. The optical spectrum of 
SDSS 0307+3848 is dominated by the secondary star and 
it presents strong emission lines. We measured the radial 
velocities from the Ha line in the SDSS-subspectra. The 
Nal and the Ha lines are in phase and the amplitudes 
were consistent with each other within the errors. 
SDSS 1006+0044, SDSS 1523+4604 and 

SDSS 1844+4120 have very cold WDs, T cS = 
7 800, 8 400, 7 600K, late secondary star spectral types 
M9, M7, M6, and are among the closest PCEBs from the 
sample, d wd = 60, 76, 73 pc respectively. SDSS 1844+4120 
is the system with the highest systemic velocity and it has 
Call H&K in absorption, i. e. the white dwarf is a DAZ. 
This photospheric absorption could be originated by some 
wind from the secondary star that is being a ccreted by 
the WD (see e.g. LTT560 (jTappert et al.ll2007D ). Looking 
at the SDSS images of SPSS 1523+4604, aka CS0749 
([Sanduleak fe Peschlll989l ). the components o f this binary 
seem to be slightly separated, also noted by iHeller et al.l 
(|2009l ). and there seems to be a third red object towards 
the North-East . We have measured an orbital period of 
9.93 hours, so it could be that SDSS 1523+4604 is a triple 
system formed by a close binary and a wide companion. 
SDSS 1143 +0009, aka WD114 0+004, is an X-ray emit- 
ting sou rce (lAgiieros et al.1 l2009t ) from the Rosat All Sky 
Survey (jVoees et al.1 1999H . The X-ray to optical flux ra- 
tio (logfx/fg) is 0.22, a value too high to be explained by 
the coronal emission of the dMe star (see iMaccacaro et al.l 
119881 : iMotch et allll998l for typical values), and the SDSS 
spectrum presents some marginal Ha emission (although 
this can be partially masked by the flux of the white dwarf). 
The high flux ratio could be explained if some wind from the 
secondary star is being accreted by the white dwarf. If so, 
some Call H&K absorption would then be expected, which 
is not visible in the spectrum. A small fraction of isolated 
white dwarfs show h ard (h i gher than 0.5 ke V) X-ray emis- 
sion (|Q'Dwver et alj 12001 IChu et al.|[200l iBih'kova et ail 
l2010l ). nevertheless the effective temperature of the white 
dwarf in SDSS 1143+0009 (~ 17 000 K) is too c old to ex- 
plain it's hard X-ray emission. Brig gs et al.1 ((2007) list other 
six PCEBs X-ray emitters: BPM 71214, RR Cae, UZ Sex, 
EG UMa, Feige 24 and WD 1541-381 . We investigated the 
X-ray fluxes of these sources within the XMM-Newton mak- 
ing use of xcat-DB^I (|Michel et al.ll2004D . None of them show 
an X-ray spectrum as hard as that of SDSS 1143+0009. It 
has a counterpart the Galaxy Evolution Explorer (GALEX) 



http:/ /xcatdb. u-strasbg.fr/ 



Nebot Gomez-Moran et al.: The orbital period distribution of SDSS PCEBs, Online Material p 2 




5 10 15 20 5 10 15 20 5 10 15 20 5 10 15 20 

Frequency [d _1 ] Frequency [d _1 ] Frequency [d _1 ] Frequency [d -1 ] 



Fig.A.l. Photometric periodogram, and photometric light curves for SDSS0853+0720, SDSS 2243+3122, 
SDSS 1611+4640 and SDSS 2112+1014 phase folded with the orbital period corresponding to the highest peak in the 
periodogram (bottom panels), P or b=0. 1503d, -P or b =0.1 1954 d, -P or b=0. 0823d and P or b=0. 0923d respectively. Double-sine 
fits are shown in red. 



database (jMartin et al.l [2005] ) . with detection only in the 
FUV flux, 17.97 ± 0.08 mag, consistent with the effective 
temperature of the white dwarf. We conclude that there is 
no straight-forward explanation of the hard X-ray emission 
of SDSS 1143+0009. 

SDSS 1352+0910 contains a hot WD, T off = 36 200 K 
and no companion visible in SDSS spectrum. Nevertheless 
emission in the Ha line was observed in the SDSS spec- 
trum, pointing to a binary nature of the system, and 
triggering further spectroscopy which revealed it to be 
a close binary. We investigated infrared surveys to learn 
more about the secondary star. We found a counterpart 
in the UKIDSS catalogue, ULAS J135228.14+091039.0. 
For details on the UKIDSS project, camera, photomet- 



ric system, calibration, pip e line 



pr ocessing and scienc e 



archive see lLawrence et all (120071) . ICasali et al.l (l2007f> 



.f (|2009h : lHamblvet all 



iHewett et all (|2006l ): iHodekin et a 

(|2008l) . Magnitudes found in the 7th data release are Y 
= 18.24 ± 0.03, J = 18.04 ± 0.03, H = +17.68 ± 0.04, 
K = 17.49 ± 0.06 mag, yielding infrared colours J-H = 
0.36 and H-K = 0.19 magnitudes. We subtracted the in- 
frared contribution of the white dwa rf and, using the em - 
pirical luminosity-mass relation from iPelfosse et all (|2000h 
we estimated a mass of 0.175 M Q for the secondary 
star. Using the spectral t ype m ass radius relation from 
iRebassa-Mansergas et all ( 20071) this is consistent with a 
spectral type of M6.5. 

SDSS 1434+5335 (SDSS J143443. 24+533521. 2) 

has a background or nearby red object 
(SDSS J143443. 19+533525.0), classified as a galaxy 
by the SDSS. The background galaxy is too faint, and 
our object is sufficiently bright so that integration time 
was relatively short and a possible contamination of our 
spectra is negligible. 

SDSS 1436+5741 and SDSS 2156-0002 present a slight 
blue excess in the SDSS spectra with respect to single 



dM stars. Spectroscopic observations carried out for these 
two systems (see Table [lj showed significant radial veloc- 
ity variation, confirming their binary nature. We explored 
GALEX with the aim of learning about the primary stars 
and found that SDSS 2156-0002 has GALEX J215614.5- 
000237 as counterpart. Only flux in the NUV was detected, 
with a magnitude of 21.99 ± 0.34, consistent with a cold 
white dwarf. 

SDSS 1559+0356 has a hot white dwarf, T cff = 48 212 K, 
and no companion is visible in the SDSS spectrum. It has a 
strong Ha emission line from which we were able to mea- 
sure the radial v elocities. SPSS 1559 +0356 was listed as a 
CV candidate bv lSzkodv et al.l (2007), but there is no obvi- 
ous sign of mass transfer in the spectrum. With its orbital 
period of 2.2 hours this system is located in the CV orbital 
period gap. If it would b e a detached CV one woul d ex- 
pect a cooler white dwarf (|Townslev fc G ansickc 2009) . We 
investigated infrared surveys to have more information for 
the secondary star, but no counterpart was found for the 
object. 

SDSS 1611+4640 has the shortest period among the 
PCEBs, with P or b = 1.97 hours, and has a cold white dwarf, 
T e ff ~ 10 300 K and a late secondary star spectral type, M5. 
SDSS 1731+6233 has a red nearby star, care was taken 
so as not to contaminate our spectra with any contribution 
from it. Unbeknownst to us at the start of our spectro- 
scopic campaign, this system had been prev iously classified 
a probable binary bv lPourbaix et all (|2005l) on the basis of 
radial velocity variations of ~ 40 kms~ . 
SDSS 2149-0717 has SDSS photometry but it was not an 
SPSS spectroscop i c targ et. It was identified as a WDMS by 
iSouthworth et al.l (|2007l) during some pilot follow-up stud- 
ies of WPMS candidates in the SPSS footprint at the UK 
Schmidt telescope with the 6dF spectrograph. The spectral 
type of the secondary star is less certain than for SPSS 
spectra since the 6dF only covers until 7 500 A. 
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SDSS 2208+1221 contains the hottest white dwarf in our 
sample, T c g ~ fOO 000 K and the earliest sp ectral type 
secondary star, a K7 star ([Silvestri et al.l [20061 ) . 
SDSS 2243+3122 has a M5 ± 1 so that it could be a de- 
tached CV (dCV) or hibernating CV, that has relaxed to 
its equilibrium radius. The SDSS spectrum presents the 
Balmer lines in emission. Radial velocities from the Ha 
emission line showed no difference in phasing from the Na I 
absorption doublet, i.e. both sets of lines must originate 
from the secondary which argues against a CV nature of 
the system. 
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Object u g r i z 

SDSS J023804.39- 000545.7 19.81 19.23 18.68 17.70 16.97 

SDSS J085336.03 + 072033.5 19.99 19.21 18.53 17.47 16.70 

SDSS J161113.13 + 464044.2 18.28 17.72 17.84 17.48 17.02 

SDSS J211205.31 + 101427.9 18.68 18.38 18.42 18.05 17.35 

SDSSJ213218.il +003158.8 18.40 18.03 18.16 17.89 17.49 

SDSS J224038.37- 093541.4 18.01 17.54 17.55 17.14 16.54 

SDSS J224307.59 + 312239.1 20.76 19.59 18.69 17.28 16.26 



Photometry: Telescope 

Ca2.2 
DuPont 
Ca2.2 
IAC80,Ca2.2 
Ca2.2 
Ca2.2 
IAC80 



duration, filter 

1.2 h, R 
3.7 h, r 
5 h, R 
19.3 h, I 
4.5 h, R 
3.5 h, R 
26.47 h, i 



R 



Table 5. Stellar and binary parameters derived for the 58 PCEBs presented in this work. 



Object 


p , 


K 

"sec 




Tvnp 


± off 


log(g) 






vd 


d w d 


Sp2 


u 2 








incl 


K ^ 


a 






Ha 






SDSS 


m 


[km s ] 


[km s J 




r t.- ] 










[M 


3) 


lP c J 


[dM] 


lP c J 




[M 


3 J 






[km s ] 


[R 
















0152 - 0058 


2.1519468(99) 


136.1(3.1) 


1.5(2.1) 


DA/M 


8773(25) 


8 


19(09) 





72(06) 


98(6) 


6 


105(49) 





20 


08) 


21 


- 22 


37(15) 


0.79(04) 





20 


09) 





88(75) 


0225 + 0054 


21.86(16) 


88.5(5.1) 


1.3(5.0) 


WD/M 


- 




- 






- 




- 


4.4 


408(147) 





32 


09) 






- 


- 







30 


10) 




- 




0238 - 0005 


5.0801(48) 


254.6(2.9) 


-13.7(2.5) 


DA/M 


20566(1537) 


7 


72 


28) 





48 


15) 


768(134) 


3 


529(104) 





38 


07) 




90 


201(71) 


1.42 


12 





39 


08) 





57 


31) 


0301 + 0502 


12.93874(47) 


176.8(2.7) 


11.1(2.1) 


DA/M 


11045(133) 


8 


41 


09) 





86 


05) 


160(11) 


4 


324(95) 





32 


09) 


08 


- 65 


65(18) 


2.82 


12) 





33 


10) 





35 


19) 


0305 - 0749 


48.4612(17) 


67.2(2.3) 


43.6(2.0) 


DA/M 


17707(410) 


7 


81 


09) 





52 


05) 


365(22) 


4 


331(97) 





32 


09) 


40 


- 46 


41(12) 


6.35 


35) 





33 


10) 





12 


07) 


0307 + 3848 


10.3266(23) 


145.1(1.7) 


-9.8(1.2) 


DA:/M 


- 




- 






- 




- 


3 


373(73) 





38 


07) 






- 


- 







39 


08) 




- 




0320 - 0638 


3.3751(18) 


287.3(3.5) 


-3.8(2.8) 


DA/M 


11173(361) 


8 


30 


25) 





79 


16) 


248(44) 


5 


423(217) 





25 


12) 


56 


- 90 


92(47) 


1.08 


10) 





26 


13) 





77 


86) 


0807 + 0724 


11.45451(17) 


119.9(5.1) 


2.7(3.0) 


DA/M 


18542(760) 


7 


98 


16) 





61 


10) 


534(56) 


3 


624(123) 





38 


07) 


41 


- 52 


74(18) 


2.51 


14) 





39 


08) 





36 


10) 


0833 + 0702 


7.344(21) 


173.9(5.8) 


88.0(4.2) 


DA/M 


15246(560) 


7 


87 


12) 





54 


07) 


386(30) 


4 


540(159) 





32 


09) 


59 


- 89 


102(31) 


1.77 


11) 





33 


10) 





44 


27) 


0853 + 0720 


3.6091(25) 


226.2(4.7) 


-10.4(3.5) 


DC/M 


- 




- 






- 




- 


3 


573(113) 





38 


07) 






- 


- 







39 


08) 




- 




0924 + 0024 


57.7( 2) 


96.0(11.0) 


1.7(6.3) 


DA/M 


19193(335) 


( 


81 


07) 





52 


03) 


335(15) 


4 


514(151) 





32 


09) 


84 


- 90 


58(17) 


7.12 


35) 





33 


10 





11 


00) 


0949 + 0322 


9.491(15) 


158.2(4.0) 


92.3(2.8) 


DA/M 


18542(737) 


•7 


79 


15) 





51 


08) 


549(52) 


4 


539(159) 





32 


09) 


60 


- 90 


98(31) 


2.07 


14) 





33 


10) 





36 


23) 


1006 + 0044 


6.7313(72) 


210.7(5.6) 


9.7(3.8) 


DA/M 


7819(61) 


8 


52 


14) 





93 


08) 


60(6) 


9 


38(1) 





12 


01) 


43 


- 48 


26(3) 


1.71 


05) 





11 


09) 





35 


32) 


1105 + 3851 


8.274(11) 


152.6(3.7) 


-9.4(2.9) 


DA/M 


10548(44) 


8 


17 


06) 





71 


04) 


206(9) 


3 


462(91) 





38 


07) 


1 ft 
41) 


- 50 


81(15) 


2.10 


07) 





39 


08) 





46 


17) 


1143 + 0009 


9.273(30) 


180.9(7.9) 


14.0(6.8) 


DA/M 


16910(349) 


7 


97 


08) 





60 


05) 


296(16) 


4 


415(122) 





32 


09) 


69 


- 90 


96(28) 


2.15 


11) 





33 


10) 





38 


22) 


1231 - 0310 


5.8493(94) 


150.9(5.3) 


36.9(4.3) 


DA/M 


10073(200) 


8 


51 


25) 





93 


15) 


162(30) 


4 


306(90) 





32 


09) 


3 ( ) 


-01 


51(16) 


1.61 


10) 





33 


10) 





63 


40) 


1300 + 1908 


7.391(13) 


138.4(4.4) 


-14.4(3.2) 


DA/M 


8673(121) 


8 


81 


19) 


1 


09 


10) 


124(21) 


4 


378(111) 





32 


09) 


28 


- 32 


40(12) 


2.09 


10) 





33 


10) 





53 


31) 


1313 + 0237 


22.31916(45) 


83.8(3.1) 


-11.0(2.3) 


DA/M 


17912(469) 


8 


09 


10) 





67 


06) 


373(25) 


4 


408(120) 





32 


09) 


.'52 


- 36 


39(11) 


3.99 


21) 





33 


10) 





22 


12) 


1316 - 0037 


9.66( 2) 


195.5(3.9) 


-12.9(3.2) 


DC/M 


- 




- 






- 




- 


3 


394(78) 





38 


07) 






- 


- 







39 


08) 




- 




1348 + 1834 


5.962( 1) 


222.0(3.8) 


-20.7(2.9) 


DA/M 


15071(167) 


i 


96 


03) 





59 


02) 


180(3) 


4 


236(69) 





32 


09) 




90 


120(34) 


1.44 


06) 





33 


10) 





56 


29) 


1352 + 0910* 


8.295(26) 


57.1(3.4) 


47.0(3.0) 


DA/M 


36154(711) 


7 


49 


11) 





45 


04) 


856(74) 




- 




- 








- 


- 






- 






- 




1411 + 1028* 


4.0160(50) 


168.3(4.3) 


42.6(4.3) 


DA/M 


30419(701) 


7 


78 


16) 





54 


08) 


829(91) 


3 


1491(294) 





38 


07) 


44 


- 57 


118(27) 


1.25 


07) 





39 


08) 





68 


29) 


1429 + 5759 


13.08486(22) 


140.5(3.9) 


-20.4(2.4) 


DA/M 


16149(109) 


8 


75 


25) 


1 


07 


13) 


329(70) 


3 


632(124) 





38 


07) 


37 


- 43 


49(11) 


3.19 


15) 





39 


08) 





38 


16) 


1434 + 5335 


104.5622(73) 


80.0(1.6) 


-40.4(1.2) 


DA/M 


21785(209) 


7 


74 


04) 





49 


02) 


153(4) 


4 


286(84) 





32 


09) 




90 


52(14) 


0.51 


46) 





33 


10) 





07 


04) 


1436 + 5741 


20.73784(71) 


120.4(2.3) 


-15.7(1.9) 


WD/M 


- 




- 






- 




- 


3 


511(101) 





38 


07) 






- 


- 







39 


08) 




- 




1437 + 5737 


20.13369(38) 


65.2(6.2) 


-2.7(4.7) 


DA/M 


17707(675) 


8 


26 


15) 





78 


09) 


319(35) 


4 


323(95) 





32 


09) 


2 1 




26(8) 


3.88 


22) 





33 


10) 





24 


14) 


1439 - 0106 


36.5426(14) 


37.7(2.9) 


9.7(2.1) 


DA/K: 


84455(5642) 


8 


11 


19) 





81 


09) 


383(65) 





826(194) 





47 


00) 


16 


- 18 


21(2) 


6.04 


15) 





49 







19 


02) 


1506 - 0120 


25.22355(39) 


134.6(4.5) 


-26.3(3.0) 


DA/M 


15601(753) 


7 


68 


19) 





45 


09) 


498(56) 


4 


526(155) 





32 


09) 




90 


95(33) 


3.99 


32) 





33 


10) 





18 


12) 


1519 + 3536 


32.81(32) 


65.2(2.3) 


-0.8(1.8) 


DA/M 


19416(246) 


7 


90 


05) 





57 


03) 


165(6) 


6 


153(71) 





20 


08) 


29 


- 30 


22(9) 


4.70 


22) 





20 


09) 





13 


11) 


1523 + 4604 


9.9317(51) 


202.5(1.9) 


-5.0(1.4) 


DA/M 


8378(49) 


8 


28 


10) 





78 


06) 


76(5) 


7 


73(40) 





15 


07) 


56 


- 63 


38(18) 


2.15 


10) 





15 


08) 





29 


30) 


1524 + 5040 


14.151167(36) 


167.5(2.3) 


-0.8(1.5) 


DA/M 


19416(292) 


8 


13 


06) 





70 


04) 


212(9) 


3 


336(66) 





38 


07) 


74 


- 90 


90(17) 


3.04 


10) 





39 


08) 





32 


12) 


1528 + 3443 


33.87(23) 


112.5(2.8) 


-79.9(2.1) 


DA/M 


26801(1022) 


7 


85 


14) 





56 


07) 


554(51) 


4 


449(132) 





32 


09) 


64 


- 90 


64(19) 


5.03 


31) 





33 


10) 





16 


10) 


1558 + 2642 


15.9016(56) 


146.7(3.4) 


-10.4(2.3) 


DA/M 


14560(367) 


8 


74 


60) 


1 


06 


31) 


283(141) 


4 


320(94) 





32 


09) 


38 


- 56 


44(17) 


3.45 


34) 





33 


10) 





32 


24) 


1559 + 0356* 


2.2662(12) 


84.3(2.7) 


-50.2(2.0) 


DA/M 


48212(2446) 


7 


99 


17) 





68 


09) 


838(106) 


_!** 


- 




- 








- 


- 






- 






- 




1608 + 0851 


9.936(25) 


172.0(4.2) 


-17.6(2.4) 


DA/M 


9844(130) 


8 


60 


15) 





98 


08) 


129(15) 


6 


201(93) 





20 


,08) 


41 


- 45 


34(14) 


2.33 


11) 





20 


09) 





35 


30) 


1611 + 0103 


7.2922(56) 


87.2(9.6) 


-11.4(4.7) 


DA/M 


10189(113) 


7 


81 


17) 





49 


10) 


230(23) 


6 


192(89) 





20 


08) 


23 


- 29 


34(15) 


1.65 


14) 





20 


09) 





34 


33) 


1611 + 4640 


1.97678(48) 


273.4(3.9) 


0.2(3.3) 


DA/M 


10307(62) 


8 


04 


09) 





63 


06) 


124(7) 


5 


274(140) 





25 


12) 


51 


- 59 


110(52) 


0.49 


03) 





26 


13) 


1 


52 


50) 


1623 + 6306 


53.5624(34) 


74.0(3.2) 


-9.8(2.1) 


DA/M 


9731(138) 


8 


64 


16) 


1 


00 


09) 


179(23) 


4 


335(99) 





32 


09) 


31 


- 34 


23(6) 


7.89 


36) 





33 


10) 





13 


08) 


1625 + 6400 


5.237716(54) 


100.5(5.2) 


-16.0(3.6) 


DA/M 


8773(76) 


8 


30 


15) 





79 


10) 


155(17) 


6 


212(98) 





20 


08) 


19 


- 22 


24(10) 


1.48 


09) 





20 


09) 





49 


44) 


1646 + 4223 


38.29(21) 


102.9(2.0) 


-14.6(1.4) 


DA/M 


17707(516) 


7 


83 


11) 





53 


06) 


399(27) 


5 


300(153) 





25 


12) 


58 


- 75 


49(23) 


5.32 


41) 





26 


13) 





13 


13) 


1705 + 2109 


19.567(61) 


94.8(3.2) 


-0.6(1.8) 


DA/M 


23886(508) 


7 


78 


07) 





52 


04) 


382(19) 


5 


276(141) 





25 


12) 


40 


- 44 


46(22) 


3.34 


22) 





26 


13) 





20 


20) 


1718 + 6101 


16.155(44) 


82.8(3.1) 


-23.7(2.4) 


DA/M 


18120(519) 


7 


83 


11) 





53 


06) 


338(24) 


4 


464(137) 





32 


09) 


32 


- 38 


49(15) 


3.07 


18) 





33 


10) 





25 


15) 


1731 + 6233 


6.4326(59) 


152.6(4.9) 


0.5(3.0) 


DA/M 


15601(488) 


7 


38 


12) 





34 


04) 


422(37) 


4 


479(141) 





32 


09) 


73 


- 90 


143(43) 


1.47 


10) 





33 


10) 





42 


25) 


1844 + 4120 


5.4166(10) 


147.2(3.1) 


-77.7(2.0) 


DA/M 


7554(6) 


7 


49 


05) 





34 


02) 


73(2) 


6 


58(27) 





20 


08) 


52 


- 57 


84(35) 


1.21 


08) 





20 


09) 





39 


33) 


2045 - 0509 


23.53152(84) 


31.1(2.7) 


-16.2(2.1) 


DA/M 


26494(1193) 


7 


69 


16) 





49 


07) 


776(83) 


3 


855(168) 





38 


07) 


14 


- 17 


24(5) 


3.94 


22) 





39 


08) 





21 


09) 


2112 + 1014 


2.2152( 1) 


- 


- 


DA/M 


19868(489) 


8 


73 


10) 


1 


06 


05) 


206(17) 


6 


199(92) 





20 


08) 






- 


0.87 


03) 





20 


09) 





97 


79) 


2114 - 0103 


9.855(14) 


136.3(2.1) 


-36.6(1.7) 


DA/M 


28064(715) 


8 


11 


12) 





70 


07) 


547(47) 


3 


591(116) 





38 


07) 


42 


- 49 


73(15) 


2.26 


10) 





39 


08) 





43 


17) 


2120 - 0058 


10.7754(26) 


57.7(2.9) 


-31.4(1.9) 


DA/M 


16336(291) 


8 


04 


06) 





64 


04) 


218(9) 


4 


316(93) 





32 


09) 


17 


- 19 


28(8) 


2.33 


,M 





33 


10) 





36 


20) 


2123 + 0024 


3.5835(74) 


149.3(6.2) 


18.3(4.6) 


DA/M 


13432(928) 


7 


31 


21) 





31 


07) 


561(79) 


6 


332(153) 





20 


08) 


42 


- 61 


94(43) 


0.85 


08) 





20 


09) 





53 


52) 


2132 + 0031 


5.3336(26) 


181.4(3.6) 


39.0(2.6) 


DA/M 


16336(303) 


7 


53 


08) 





39 


03) 


343(18) 


4 


530(156) 





32 


09) 




90 


148(43) 


1.33 


07) 





33 


10) 





50 


28) 


2149 - 0717 


15.453(11) 


100.7(0.9) 


-1.1(0.6) 


DA/M 


14393(113) 


8 


11 


03) 





68 


02) 


144(3) 


2 


220(52) 





44 


!n) 


39 


- 41 


64(1) 


3.21 


02) 





44 


10) 





32 


01) 


2156 - 0002 


12.07015(25) 


94.7(2.6) 


-39.2(1.5) 


DA:/M 


















6 


177( 82) 





20 


08) 















20 


99) 








2208 + 1221 


45.67(15) 


69.5(1.4) 


-25.7(1.1) 


DA/K 


99575(5434) 


9 


02 


34) 


1 


24 


13) 


383(152) 


-1** 






























2216 + 0102 


5.0487(51) 


193.6(4.0) 


-11.3(2.8) 


DA/M 


12536(1541) 


7 


62 


28) 





41 


14) 


229(37) 


5 


360(184) 





25 


12) 


58 - 


120(70) 


1.30 


18) 





26 


13) 





46 


59) 


2240 - 0935 


6.2541(27) 


160.7(1.9) 


-15.7(1.4) 


DA/M 


12536(1036) 


7 


62 


14) 





41 


06) 


210(18) 


5 


231(118) 





25 


12) 


55 


- 83 


99(49) 


1.47 


14) 





26 


13 





41 


44) 


2243 + 3122 


2.8703(59) 


184.6(5.6) 


-19.4(4.3) 


DC/M 


















5 


171(87) 





25 


12) 













26 


13) 






2311 + 2202 


13.9224(17) 


114.8(2.7) 


-1.3(2.0) 


DA/M 


10189(110) 


9 


06 


09) 


1 


23 


05) 


115(12) 


3 


474(93) 





38 


07) 


29 


- 30 


35(6) 


3.30 


08) 





39 


08) 





39 


14) 


2318 - 0935 


60.8094(27) 


98.2(3.1) 


5.0(1.9) 


DA/M 


22550(603) 


7 


75 


09) 





50 


05) 


460(29) 


3 


791(156) 





38 


07) 




90 


74(15) 


7.49 


33) 





39 


08) 





11 


04) 



Notes. Radial velocities are measured from the Nal doublet, unless marked with * in which case it comes from the Ha emission line. K w d are the calculated values in 
Sect. I4.2~3l Errors in the last digits are given in parenthesis. ** No spectral type could be derived from fitting the spectrum to template spectra. 
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Fig. 3. SDSS spectra of the 58 systems with orbital period measurement. Systems are sorted in right ascension from left 
to right and from top to bottom. 
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Fig. 4. Phase folded radial velocities curves and sine fits to the data (dashed lines). As for Fig. [3] systems are sorted in 
right ascension from left to right and from top to bottom. 
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Fig. 5. SDSS spectra of the 58 systems with orbital period measurement. Systems are sorted in right ascension from left 
to right and from top to bottom. 
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Fig. 6. Phase folded radial velocities curves and sine fits to the data (dashed lines). As for Figure [3] systems are sorted 
in right ascension from left to right and from top to bottom. 



